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ABSTRACT 

We study accretion in binaries hosting an intermediate mass black hole (IMBH) of 
^ . ~ 1000 Mq, and a donor star more massive than 15 M,q. These systems experience an 

' active X-ray phase characterized by luminosities varying over a wide interval, from 

£SJ | < 10 36 ergs _1 up to a few 10 40 ergs _1 typical of the ultra luminous X-ray sources 

■ (ULXs). Roche lobe overflow on the zero-age main sequence and donor masses above 

r*"** ' 20 M can maintain a long-lived accretion phase at the level required to feed a ULX 

source. In wide systems, wind transfer rates are magnified by the focusing action of 
j^j ■ the IMBH yielding wind luminosities Z 10 38 ergs -1 . These high mass-IMBH binaries 

can be identified as progenitors of IMBH-radio pulsar (PSR) binaries. We find that the 
formation of an IMBH-PSR binary does not necessarely require the transit through a 
<"*> |' ULX phase, but that a ULX can highlight a system that will evolve into an IMBH- 

PSR, if the mass of the donor star is constrained to lie within 15 to 30 M Q . We show 
that binary evolution delivers the pre-exploding helium core in an orbit such that 
after explosion, the neutron star has a very high probability to remain bound to the 
IMBH, at distances of 1-10 AU. The detection of an IMBH-PSR binary in the Milky 
Way has suffered, so far, from the same small number of statistics limit affecting the 
population of ULXs in our Galaxy. Ongoing deeper surveys or next generation radio 
telescopes like SKA will have an improved chance to unveil such intriguing systems. 
jJJ ■ Timing analysis of a pulsar orbiting around an IMBH would weigh the black hole in 

the still uncharted interval of mass around 1000 M Q . 

Key words: ULX — galaxies: neutron star — IMBH — X-rays: binaries — X-rays: 
galaxies 
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1 INTRODUCTION 

Recent high resolution X-ray imaging and spectroscopic 
studies with Chandra and XMM, have led to the discov- 
ery of a large sample of a new class of compact sources 
with luminosities in the interval between 3 x 10 39 erg s _1 
and 10 41 ergs -1 , that are in excess of the Eddington limit 
of a stellar-mass black hole of 20 Mq (Fabbiano 1989; see 
Mushotzky 2004 for a critical review). These sources can find 
a simple interpretation in the hypothesis that intermediate 
mass black holes (IMBHs) exist with mass 10 2 M - 10 4 M o 
accreting from a companion star in binary systems (Fab- 
biano 1989; Mushotzky 2004; Miller & Colbert 2004). The 
detection of a cool-disc thermal spectral component in a 
number of ULXs (Miller et al. 2003; Zampieri et al. 2004; 
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Miller, Fabian & Miller 2004; Kaaret et al. 2004; Cropper 
et al. 2004; Dewangan et al. 2004), the properties of the 
optical and radio counterparts (Koerding, Colbert & Fal- 
cke 2005; Liu et al. 2004; Zampieri et al. 2004; Kaaret et 
al. 2004; Soria et al. 2005; Miller, Mushotzky & Neff 2005), 
and the timing behaviour of at least one source in the star- 
burst galaxy M82 (Strohmayer & Mushotzky 2003; Fiorito 
& Titarchuk 2004), support this view. The IMBH hypothe- 
sis however does not represent the only possibility to explain 
the emission of ULXs, since mechanical beaming working in 
a thick disc around a conventional stellar-mass black hole, 
or Doppler boosting from a jet in a microblazar could pro- 
duce the same range of observed luminosities (King et al. 
2001; Koerding et al. 2002; Mushotzky 2004). Moreover, in 
a recent work of Rappaport, Podsiadlowski & Phfal (2005) 
the authors use binary evolution calculations to show how 
the largest part of the ULX population may be explained 
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with a stellar mass black hole emitting at a super Edding- 
ton rate of ~ 10 without the requirement of an IMBH (see 
also Podsiadlowski, Rappaport & Han 2003 and Phfal, Pod- 
siadlowski & Rappaport 2005 for an extended study on the 
evolution of stellar mass black hole binaries). This is also 
in agreement with another study of King & Dehnen (2005) 
where the authors claim that the luminosities of a large sam- 
ple of ULXs could be explained using helium enriched mat- 
ter and mechanical beaming with a stellar mass black hole. 
But, all the alternatives to the IMBH hypothesis meet with 
strong difficulties when the luminosity of a ULX is in ex- 
cess of ~ 10 40 ergs _1 , since, in this case, beaming under 
rather extreme conditions should be at work to match with 
the observations, or super Eddington factors greater than 
~ 10 could be difficult to achieve. On the other hand, al- 
though the hypothesis of an IMBH explains naturally many 
of the observational clues, it clashes with the problem of 
providing a viable mechanism of formation. Until now, two 
possibilities have been proposed: the formation of an IMBH 
through runaway collisions among massive stars undergoing 
fast dynamical segregation, in the core of a dense super star 
cluster (see Porteigies Zwart et al. 2004; Gurkan, Freitag & 
Rasio 2004), or the wandering of an IMBH, relic of a zero 
metallicity population III star (Abel et al. 2000). In the first 
case, the giant star that forms in the core of the dense star 
cluster, collapses into an IMBH. This can occur in a star- 
forming region, and the result is based on very large detailed 
N-body simulations (Portegies Zwart et al. 2004a). The case 
of wandering IMBHs relic of the early assembly of haloes in 
a currently starforming galaxy is uncertain, in particular the 
capture of gas or of a star to ignite accretion (Volonteri & 
Perna 2005). 

In order to solve the controversy on the real existence of 
IMBHs in ULXs, the only secure route would be the determi- 
nation of the optical mass function, similar to the procedure 
for the stellar-mass black holes in the Milky Way (Orosz et 
al. 2004). This is difficult however, since ULXs are distant 
sources hosted in external (starburst) galaxies for which the 
optical identification of the companion is troublesome, and 
even in the lucky circumstance of good identification (see 
Zampieri et al. 2004; Kaaret et al. 2004; Soria et al. 2004; 
Liu et al. 2004; Miller et al. 2004) the optical spectrum is 
too noisy to allow a mass estimate of the black hole. 

In this paper we propose an alternative way to discover 
and weigh a IMBH: it uses the detection of a young ra- 
dio pulsar around an IMBH. If ULXs are indeed accreting 
IMBHs in binaries, their evolution would end with a radio 
pulsar around the IMBH, if the mass of the donor star is in 
the range which avoids the formation of a stellar-mass black 
hole or of a white dwarf. The detection of a radio pulsar 
orbiting around an IMBH would provide, through timing, 
an unambiguous measure of its mass. 

At present, no young radio pulsars have been seen orbit- 
ing a heavy invisible companion, but in the near future the 
search of these hypothetical IMBH-pulsar systems (IMBH- 
PSR hereon) can be extended to nearby galaxies, with LO- 
FAR (Rottgering 2003) and SKA (Cordes et al. 2004). An 
IMBH-PSR system would lead to the discovery of black 
holes in the still uninvestigated interval of masses between 
100 M Q - 1O 4 M . Stellar -mass black holes around pulsars 
have long been considered to be the "holy grail" of com- 
pact star binaries. The systems studied here are even more 



intriguing, given their importance in discovering a new un- 
explored mass range crucial for cosmology (Madau & Rees 
2001). The discovery of a pulsar orbiting around a stellar- 
mass black hole is expected in the next years, on the basis of 
theoretical considerations concerning their population, and 
on the observability of the radio pulsar signal (Lipunov et 
al. 1994, Sigurdsson 2003, Lipunov et al. 2005). On similar 
lines, there is also the hope of detecting, despite the large 
interstellar electron densities, a radio pulsar orbiting around 
the massive black hole of ~ 3 x 10 6 M Q hosted in the core of 
our Galaxy (Pfahl & Loeb 2004). Thus, the issue of pulsars 
around black holes is becoming of paramount importance. 

In this paper we will assume the hypothesis of the for- 
mation of IMBHs in young dense star clusters, and will start 
our evolution study just after the formation/capture of a 
high mass star around the IMBH, which could be the pro- 
genitor of an IMBH-PSR system. In this framework, Hop- 
man et al. (2004) considered the possibility that a passing 
star is tidally captured by the IMBH in a stable, close, not 
plunging orbit. Mass transfer may initiate, after circulariza- 
tion, while the star is on the main sequence or evolving away 
from it. Dynamical capture of a massive star by an IMBH is 
a further possibility, as shown by Baumgardt et al. (2004). 
In an exchange interaction of a binary star, the IMBH can 
acquire a companion, likely a massive star, given that the 
IMBH forms in a mass-segregated environment, where stel- 
lar encounters play an important role. 

The observational appearance of binaries hosting an 
IMBH has only been partly explored, and mainly in the con- 
text of ULXs. Portegies Zwart et al. (2004b) studied the evo- 
lution of an IMBH of 1000 M accreting from a donor star 
of mass between 5 to 15 M (see also Kalogera et al. 2004 for 
another binary evolutionary study). Typically, light donors 
( ;£ 5 M ) do not transfer mass at a sufficiently sustained 
rate to produce sources as bright as ULXs; only the high 
mass stars (^1OM ) on the main sequence or beyond can 
provide luminosities in excess of 10 40 ergs -1 . Thus, IMBHs 
cannot easily be identified on the basis of their X-ray ac- 
tivity: they can display a rather wide range of luminosities, 
depending on the mass of the donor star, and on the mass 
transfer mechanism, as it will be shown in this paper. Only 
when they outshine as bright ULXs can they become visible 
over the stellar-mass black holes. 

In this paper we address a number of issues, in the hy- 
pothesis that an X-ray and/or a ULX phase preceeds the 
one in which the system appears as an IMBH-PSR binary: 

• What are the characteristics of the X-ray phase and of 
mass transfer when considering massive donors around an 
IMBH? 

• Is the formation of a radio pulsar likely? 

• What is the probability that an asymmetric kick, im- 
parted to the neutron star at birth, will unbind the system? 

• On which orbits the radio pulsar is released after the 
supernova explosion? Can Doppler effects hamper the de- 
tection of the radio signal? 

In order to address these questions, we explore in §2 wind fed 
accretion (WFA), and accretion through Roche lobe over- 
flow (ROLF), when considering donor stars with masses in 
excess of 15 M . We infer the typical luminosities of these 
high-mass IMBH binaries and the fate of the mass tranfer- 
ring star. In §3 we explore natal kicks in such exotic massive 
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binaries, and determine probabilities of survival of the bi- 
nary. Taking into account the effect of the natal kick, wc 
also compute the typical eccentricity and semimajor axis of 
the neutron star orbit, after the supernova explosion in the 
binary. In §4 we simulate the observability of an IMBH-PSR 
in the Galaxy, considering the parameters of the recent sur- 
vey for radiopulsars in the Galactic disc using the Parkes 
radio telescope. In §5 we discuss our results and outline our 
conclusions. 



2 ACCRETION ONTO IMBHS 

Here we explore the evolution of binaries with an IMBH of 
1000 Mq and stars having mass M heavier than 15 M©, 
so that we can be confident that, despite mass loss, the 
donor star does not evolve into a white dwarf. The he- 
lium core mass necessary to form a neutron star is in the 
range 2.8M© < Mh c < 8M©, although some uncertainites 
exist in these upper and lower limits (Tauris & van den 
Heuvel 2004). In each evolutionary sequence, obtained by 
varying the mass and the initial orbital separation, the star 
is evolved until carbon ignition. The mass transfer episodes 
are either driven by RLOF during the main sequence phase; 
by RLOF during the rapid expansion of the star when as- 
cending the giant branch; or by WFA when the system is 
detached. This third possibility has always been overlooked 
in previous studies, though it seems relevant to understand 
the different phases of an accreting IMBH. 

The evolution equations for the binary systems are 
solved using an updated version of the Eggleton code (Eggle- 
ton 1971; Pols et al. 1995) with a modified mass transfer 
rate for the WFA case, as explained later in the text. All 
the stars considered have Population I chemical composi- 
tion (Y=0.28, Z=0.02), radii R and wind losses M w obtained 
according to de Jager et al. (1988). The mixing length pa- 
rameter is set to a — 2.0, and the overshooting constant 
8 v = 1-2 (Pols et al. 1998). Using general arguments, we 
then focus attention on the remnant left in the evolution 
to study the detectability of a young neutron star revolving 
around an IMBH. 



2.1 Wind Fed Accretion 

Consider a donor main sequence star of mass M above 
15 M© , moving on a circular orbit around an IMBH of mass 
Mbh. Let a be the relative separation, 
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the stellar wind velocity, set equal to the escape velocity 
from the stellar surface. 

WFA dominates when the star is underfilling its Roche 
lobe radius, i.e., when the binary separation a exceeds orl, 
denned as 



R[0.6q 2/ '- +\n(l + q 1/s )] 
0.49g 2 / 3 



(3) 



where q = M/M B u- 

In high mass X-ray binaries, the mass of the donor usu- 
ally exceeds the mass of the black hole, hence, the wind 
particles mainly accrete onto the black hole from unbound 
orbits, since V OT h <§C V w . Capture occurs within a cylinder of 
radius comparable to obh ~ 2GM B h/V3 « (Mbh/M) R. In 
the fluid approximation, this leads to an accretion rate onto 
the hole of 



M BH ~ 1/4 (a BH /a) 2 (Kd/Kv)M» 



(4) 



(Shapiro & Lightman 1976, Shapiro & Teukolsky 1983), 
where V le \ is the wind speed relative to the black hole 



cos a, 



(5) 



with a £ [0, 7r] the angle between the directions of V or b and 
Vw Equation (4) is valid when a 3> (Zbh, and Vorb <C Vw, so 
that Kei ~ Vw. 

In a binary with an IMBH, we may have obh Z a (or 
equivalently Vorb^Vw), i.e., the particles released from the 
donor star surface can have kinetic energies such to be bound 
to the IMBH. This effect may lead to an accretion rate Mbh 
comparable to the wind mass loss rate M w • In particular, if 
all the wind particles are deep in the potential well of the 
IMBH, their total energy K + U = V^/2 - GM BH /a is 
negative. This implies a full wind gravitational focussing by 
the IMBH, and an X-ray luminosity at its highest value. The 
limiting orbital separation for this to occur corresponds to 
the case K + U = for the most energetic wind particles, 
i.e., the ones having a = 0. This yields V or b = [1 + v2]Vw, 
and accordingly, a binary separation 



GMbh 



(1 + V2) 2 Vi 
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Below a m in the wind is entirely accreted. When a — n, the 
condition K + U = corresponds to those particles having 
the minimum relative speed to remain bound to the IMBH. 
This defines a critical orbital velocity V or b = [a/2 — 1]Vw and 
in turn a limiting orbital separation 



GMbh 



{y/2-l)*V* 



M 



(7) 



above which all wind particles are unbound. In the interval 
a m in < a < dmax wind particles with K + U > escape the 
gravitational field of the IMBH and this defines a loss cone 
of solid angle Q — 2ir(l — cos a), with a = cos _1 [(V 2 b — 
K 2 )/(2VwVorb)]. Under these conditions, the accretion rate 
onto the IMBH can be approximated as 



Mbh ~ Max 



4tt 



) M w ;M!h 



(8) 



For the mass transfer rate of eq. (8) even a very high mass 
star can provide the ULX luminosities during the main se- 
quence, without the requirement that the donor is cross- 
ing the supergiant phase. In Figure 1 we illustrate this ef- 
fect, plotting accretion rates Mbh and Mbh f° r an IMBH 
of 1000 M© and a donor star of 30 M© during the main se- 
quence. The probability to observe this system is strongly 
increased due to the larger time spent by the star on the 
main sequence. 
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Figure 1. Mbh (solid line) and Mjjg (dashed line) versus a for 
a donor of 30 Mq and an IMBH of 1000 Mq. The mass transfer 
rates are in units of the wind mass loss rate. The accretion rate 
onto the IMBH is the highest value between the two. The figure 
is explained in the text. 

In Figure 2, we show a m i n and a max as a function of 
the mass of the IMBH, compared to anh, i.e., the distance 
at which the donor star overfills its Roche lobe. The light 
gray area corresponds to orbital separations where WFA 
is in the Shapiro-Lightman limit, while in the middle gray 
area, the wind is enhanced by the action of gravitational 
focussing by the IMBH, and in the heavy gray area, it is 
fully captured, leading to a conservative system. The black 
area is the RLOF zone. 

In our WFA model, the possible formation of a cir- 
cumbinary disc around the donor star might change the 
quantity of wind accreted, although the order of magnitude 
remains unchanged (see van den Heuvel 1994). We have 
followed the evolution of the binary under WFA, starting 
from a separation of 1 AU, imposing angular momentum 
and mass loss from the fraction of the wind escaping from 
the binary. The star is non rotating and the wind leaves the 
donor isotropically, carrying away, in the centre mass refer- 
ence frame, a specific angular momentum J = h(M w — Mbh) 
with h the specific angular momentum of the orbit (see 
Soberman, Phinney & Van den Heuvel 1997 for a detailed 
discussion). We have computed the accretion rate onto the 
IMBH using eq. (8) and a WFA luminosity adopting an ef- 
ficiency factor of ~ 0.1 characteristic of disc accretion, since 
the captured wind has high enough angular momentum to 
form a disc. 

The run of the WFA luminosity as a function of time 
is shown in Figure 3 for an IMBH of 1000 Mq in a binary 
with a 20-30-50 Mq star, respectively. The luminosity is en- 
hanced relative to the value that would be predicted from 
a stellar-mass black hole under equivalent conditions. With 
time the WFA luminosity rises mainly because of the instrin- 
sic increase of the wind loss rate when the star is transiting 
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Figure 2. (Irl (dashed line), a m i n (lower solid line) and a max 
(upper solid line) versus Mbh for a donor of 30 Mq on the ZAMS. 
WFA prevails above the dashed line. If the donor star is inside the 
gray shaded area, WFA is enhanced by the gravity of the IMBH. 
Black area (flmin ^ ^rl) corresponds to RLOF, while the heavy 
gray area is for accretion of the entire wind lost by the donor star. 
Mild gray refers to a partial focussing of the wind, while the light 
gray is the Shapiro-Lightman regime. The condition a m - ln = orl 
defines a minimum IMBH mass for which WFA occurs under full 
capture of the stellar wind, giving maximum WFA luminosities. 
For a 30 Mq donor this minimum IMBH mass corresponds to 
~ 3000 M . 

from the zero age to the terminal age main sequence. Only 
very massive stars can provide mean luminosities in excess 
of 3 x 10 39 ergs -1 , the threshold defining ULX sources. 

2.2 Accretion through Roche lobe overflow 

Accretion via RLOF occurs during the main sequence when 
the star is close enough to the IMBH to fill its Jacobi surface, 
or during the rapid expansion of the star when ascending the 
giant or supergiant branch. In this section we explore the 
evolution for the former case, i.e., for separations ranging 
between 0.18-0.25 AU for a 1000 Mq IMBH, and 15-50 M 
companion stars. 

Roche lobe contact remains stable along the entire 
evolution and gives luminosities above 10 39 ' 5 ergs _1 when 
the donor star has a mass M>15M . Luminosities of 
10 40 ergs _1 , typical of the brightest ULXs, need a star at 
least as massive as 20 Mq , as illustrated in Figure 4. Adopt- 
ing the Dubus criterion for assessing the stability of the 
disc against the thermal ionization instability (Dubus et al. 
1999), we find that RLOF onto an IMBH is stable, under 
our conditions. This ensures that high mass IMBH binaries 
can remain bright for a relativley long time, between a few 
million to ten million years. This is opposite to the case 
of IMBH binaries with donor stars having masses < 10 M 
as pointed out by Portegies Zwart et al. (2004b), for which 
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Figure 3. From the top to the bottom lines: WFA luminosity 
against time for donors with masses of 50 Mq, 30 Mq, and 20 Mq 
respectively, and a black hole's mass of 1000 Mq. The initial or- 
bital separation is 1 AU. The dashed line is the limit to obtain an 
ULX. The mass transfer begin on the ZAMS and the integration 
is stopped when the star reaches the giant branch. At this point a 
RLOF phase starts, producing the sharp spikes in the plot, with 
the exception of the donor of 50Mq for which the code is stopped 
at the end of the main sequence. 



the luminosities are transient and the X-ray phases last for 
a longer time, comparable with the donor lifetime. Within 
the limits of the available time span coverage of ULX ob- 
servations, the apparent persistency in their emission may 
be an indication that donor stars have masses in the range 
considered in this paper. 

The high mass transfer rates that establish during 
RLOF dramatically change the mass of the star. Figure 5 
shows the donor mass as a function of time. At the end of the 
stellar evolution, the core mass is almost equal to the mass 
of the whole star, since the mass transfer strips almost com- 
pletely the entire stellar envelope during the RLOF phase. 
Donors with mass between 15 Mq and ~ 25Mq on the zero 
age main sequence likely end their lives with a core mass be- 
tween 3.3 and 7M , probably producing a neutron star. The 
case of a 30Mq could be considered the mass limit between 
the neutron star and the black hole formation, since its final 
core mass is ~ 9 Mq . The fate of the heaviest star here con- 
sidered (M>30Mq) is to collapse into a stellar-mass black 
hole. 

If RLOF is the mechanism leading to a ULX during the 
main sequence, the donor star is expected to lie at an initial 
distance of a S 1 AU from the IMBH. After leaving the main 
sequence, the star ascends the giant branch, while keeping 
on overflowing its Roche lobe. This generate the spikes in 
Figure 4 giving very high luminosities. The typical values 
of a reached at the end of the evolution are in the range of 
~ 2 - 5 AU. 

If the initial orbital separation of a donor star is a > 1 



Figure 4. Luminosity against time, for donors with masses of 
15,20, 25, 30, 50M Q , respectively, and for an IMBH of IOOOMq. 
The mass transfer begins near the zero age main sequence and 
the integration is stopped before carbon ignition. Stellar masses 
in excess of ~ 20Mq can produce luminosities above 10 40 ergs -1 
during the main sequence. All binaries appear to be persistent 
sources. The dashed line represents the ULX limit. Note that 
there is only one spike for each donor, corresponding to the giant 
phase. The envelope of the donor is almost completely depleted 
during this phase and the rc-cxpansion after the helium burning 
is too small to produce another contact phase. 



AU, RLOF is avoided on the main sequence phase. The star, 
in exiting the terminal age main sequence, may later over- 
fill its Roche lobe due to the sizable increase of its radius 
when ascending the giant branch, initiating a phase of ac- 
cretion. Mass transfer and angular momentum conservation 
will drive the star further away from the IMBH, up to typical 
distances of 8 — 10 AU. 



3 NEUTRON STAR KICKS AND THE 

SURVIVAL OF THE IMBH-PSR BINARY 

In the RLOF case, the continuous increase of a, implied by 
conservative mass transfer, would bring the donor to dis- 
tances of several astronomical units. We may ask whether 
the explosion and the asymmetric kick that accompanies the 
formation of the neutron star can unbind the system. Given 
the high inertia of the IMBH, a neutron star has a greater 
chance to remain bound than in the case of a lighter black 
hole companion. 

To keep a binary bound, we need a kick velocity below 
the limit given by the equation (Brandt & Podsiadlowski 
1994, Willems, Kalogera & Henninger 2004) 



G (M BH + M) 



ao 



1/2 



/ M BH + M NS \ 
+ V Af BH + M J 



1/2' 



(9) 



where Mns is the mass of the neutron star formed after 
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Figure 5. Donor mass versus time during the whole life of the 
star for RLOF systems. The mass loss is very rapid for a high mass 
star and decrease mildly for a lighter star. The plot shows binaries 
with donor star masses of 15, 20, 25, 30 and 5OM0 orbiting around 
a IOOOMq IMBH. The very rapid short decrease near the end of 
the curves is due to RLOF on the giant phase. Note that the initial 
mass corresponds to an age of 1 Myr as we are using logarithmic 
scale on the time axis. 



the supernova explosion, and M and ao are the mass of the 
donor star and binary separation, just prior the supernova 
explosion. In the presence of the IMBH, the total mass re- 
mains almost constant before and after the explosion, and 
so the most important parameter is the orbital separation 
ao. To obtain a conservative estimate, we use our largest 
value of ao which is about 10 AU. The kick velocity limit is 
around 700 km s -1 . 

To calculate the probability of survival of the IMBH- 
PSR system, we adopt an analytical expression for the dis- 
tribution of the kick magnitudes (Hansen & Phinney 1997, 
HP distribution) 



p(Vfc) = 



-Vila' 



(10) 



where a is the velocity dispersion of each of the components 
of the kick velocity. 

Although there is not full agreement on the form of the 
kick magnitude distribution, it is rather well established that 
speeds larger than ~ 700 km s" 1 are attained only in extreme 
rare cases (see also Lyne & Lorimer 1994, Arzoumanian et 
al. 2002, Hobbs et al. 2005). Therefore we are confident that 
our calculations are not strongly affected by the particular 
choice of the distribution, since what is interesting here is 
the existence of a cutoff at high velocities more than the 
specific form of the distribution itself. 

We use the following equation to calculate the proba- 
bility of survival of our system, obtained by integrating the 
HP distribution between zero and our kick velocity limit: 



-oi Q. 




Figure 6. Distribution of post-supcrnova systems over a and e, 
normalized as in equation (14) of Kalogcra (1996), calculated for 
P ~ 1 and £ = 0.3. 
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Using a value for 14, i of 700 kms~ and a value for a of 152 
kms -1 (Hobbs et al. 2005), we infer a probability of survival 
greater than 99%. Therefore after the supernova explosion 
we can be confident that the system is still bound. 

After the explosion the orbital parameters change. To 
obtain the distribution of post-supernova orbital separations 
and eccentricities, we introduce, following Kalogera (1996), 
the three dimensionless parameters that are necessary to 
constrain the properties of the binary after the explosion: 



P = 

£ = 



a//ao 

Mbh + Mns 



(12) 
(13) 
(14) 



M BH + M 
a 

Vorb 

where a is the pre-explosion binary separation and aj its 
post-explosion value, f3 ~ 1 the binary total mass ratio, after 
and prior explosion. 

Typical pulsar ID speeds peak at 152km s _1 (Hobbs 
et al. 2005). Given the high relative orbital velocity (eq. 1) 
with respect to a, £ < 1 in our case. Figure 6 shows the prob- 
ability distribution over a and post-encounter eccentricity e, 
computed for (3 ~ 1 and £ = 0.3, while Figure 7, and Fig- 
ure 8 give the distribution over e (integrated over a) and a 
(integrated over e) respectively, for selected values of £. 

Post- supernova binaries populate only a restricted area 
of the a—e plane giving minimum acceptable post- supernova 
separations a/ = (l/2)a for £ ~ 1 and a/ = ao for £ << 1 
according to the fact that the post-supernova orbit must 
include the position of the two stars just prior to explosion. 
Because we have found in Section 2.2 that the likely pre- 
supernova orbital separation falls in the range 2-10 AU, final 
values of a/ range between ~ 1 — 10 AU. 

The small eccentricity of these systems is a direct conse- 
quence of the very high mass of the IMBH, which produces 
higher orbital velocities than in normal binaries, and a value 
of £ < 1 also for the largest orbital separations. Another 
consequence of the large mass of the black hole is that the 
binary can be always defined as an hard one (Heggie 1975) 
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Figure 7. Distribution of post-SN systems over eccentricites, cal- 
culated for different values of £ and for fi ~ 1. For £ = 1 the sys- 
tems are highly eccentric, while they became less eccentric with 
decreasing £. 



and thus the dynamical interactions in the star cluster tend 
to decrease the orbital separation of the system, although 
there are some uncertainities in the real effect of this pro- 
cess as stated by Hopman & Portegies Zwart (2005). 

Neglecting this uncertain effect, the production of grav- 
itational waves seems to be unimportant, because the pulsar 
is revolving on a mildly eccentric orbit at a large separation. 
The in-spiral time for emission of gravitational waves is ex- 
tremely large: £gw £ 3 x 10 10 yrs for e = 0.4, and only for 
very high eccentric systems (e = 0.9) it is £gw ~ 6 x 10 7 
yrs. 

The lifetime of a super star cluster is about 10 s yrs (see 
Hopman & Portegies Zwart 2005), which is comparable to 
the lifetime of the binary even if the eccentricity produced 
by dynamical interactions becomes ~ 0.9. Thus we have a 
very small probability of a merger in a super star cluster by 
gravitational waves emission, and the observability in the 
radio band of our system is not compromised by this effect. 
This is in agreement with the results of Hopman & Portegies 
Zwart (2005) which found that with an IMBH of ~ 1000 M , 
only a fraction around 2% would merge in an Hubble time. 



4 SIMULATING PULSARS AROUND IMBH 

The search of pulsars in binary systems is more difficult than 
for isolated pulsars, because of changes in the line-of-sight 
velocity due to orbital motion cause varying observed pulse 
periods. Most of the algorithms for searching for periodic- 
ities in a time series use a Fast Fourier Transform (FFT) 
to produce a frequency spectrum; a varying received pul- 
sation frequency spreads the detection over many spectral 
bins. This reduces the signal-to-noise ratio (s/n) of a detec- 



Figure 8. Distribution of post-SN systems over dimcnsionlcss 
orbital separation, calculated for different values of £ and for (3 ~ 
1. For £ = 1, the most probable post-SN semimajor axis is about 
0.5ao, while for £ << 1 the most probable semimajor axis is 
similar to the prc-SN orbital separation. 

tion with respect to that of an isolated pulsar with the same 
intrinsic flux and rotational period. 

We have investigated the consequencies of this effect on 
the probability of detecting a pulsar orbiting an IMBH dur- 
ing the Parkes Pulsar Multibeam survey (Manchester et al. 
2001), i.e. the deepest and most successful large scale survey 
for pulsars in the Galactic disc completed so far, resulting in 
the discovery of more than 700 new sources (Faulkner et al. 
2004) . We have generated many sets of simulated time series 
containing a pulsating signal (with a duty cycle of 5%, as 
statistically appropriate for non-recycled pulsars) of a pulsar 
of 1.4 M orbiting an IMBH of 1000 M Q . The other binary 
parameters have been chosen to span the most likely values 
of the distribution plotted in Figure 6, i.e. orbital separations 
in the range 1-10 AU (corresponding to orbital periods from 
~ 10 days to ~ 1 yr) and eccentricity from to 1 in steps of 
0.05. The sampling time (0.25 ms) and the total number of 
samples (2 23 ) are the same as the observations of the Parkes 
Pulsar Multibeam survey. Spin periods have been chosen to 
be 1000 ms, 100 ms and 10 ms (the latter being a rather 
extreme case for typical non-recycled pulsars, which would 
be expected to be hosted in a binary with an IMBH). For 
each of the parameters, we have produced eight time series, 
spanning four different values of the longitude of the perias- 
tron of the elliptical orbit (0, 90, 180 and 270 deg) and two 
values of the orbital phase (quadrature or conjunction) for 
simulated observation. 

After creating this large set of time series, we first ap- 
plied the standard search (which maximizes the sensitivity 
for detecting isolated pulsars) and, later, the "stack search" , 
which involves splitting each time series into 16 segments, 
performing an FFT on each of them. Then the resulting 
spectra are summed applying a constant shift between fre- 
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quency bins of adiacent spectra, the shift corresponding to 
the change in the apparent spin frequency between one seg- 
ment of observation and the next (Faulkner et al. 2004). 
Many different bin shifts are explored in order to find the 
ones which maximize the spectral signal-to-noise of a given 
periodicity in the whole time series. At the sacrifice of a cer- 
tain loss of sensitivity (due to the incoherent sum of spec- 
tra) this algorithm considerably improves the possibility of 
detecting Doppler distorted signals at a reasonable compu- 
tational cost. 

For each simulated time series, we finally compared the 
signal-to-noise ratio (s/n) obs at which the simulated pulsa- 
tion period has been recovered by the search analysis with 
the signal-to-noise ratio (s/n) so i of the same signal if it were 
been emitted by a solitary pulsar. For spin period P > 100 
ms, it results that rj = (s/n) obs /(s/n) so i > 0.8 for all cases, 
but those involving e > 0.7 and the closest orbital separa- 
tions. The application of the stack search allows to mantain 
a relatively good sensitivity (rj £ 0.6) at most orbital phases 
and eccentricities even for a 10 ms pulsars. In summary, only 
the combination of very extreme eccentricities e > 0.9, un- 
favorable orbital phases (longitude of periastron ~ 180 deg 
and an observation performed at about the epoch of perias- 
tron) and orbital separations in the lower end of the assumed 
range decrease the sensitivity of the search below rj = 0.5. 

For a disc structure with thickness much smaller than 
its radial extension, rj roughly represents the ratio between 
the galactic volumes explored to search for the two types of 
sources; thus, assuming a disc-like and uniform distribution 
for both pulsars, i.e., those in IMBH-PSR systems and the 
isolated pulsars, rj indicates approximately the relative de- 
tection efficiency. In summary, given the orbital parameters 
of the putative pulsars orbiting IMBHs, under the conditions 
explored, we conclude that their detectability by the Parkes 
Multibeam Pulsar survey has only been modestly affected 
by their orbital motion. Other biases related to large-scale 
pulsar surveys (discussed in §5.2) play the major role in 
selecting the observed population from the intrinsic popula- 
tion of IMBH-PSR systems. 



5 DISCUSSION 

5.1 The formation of an IMBH-PSR and its 
X-ray luminosity 

The formation of an IMBH-PSR binary system relies mainly 
on four working hypothesis: that an IMBH (i) forms in a 
young core-collapsed dense star cluster, (ii) acquires, prefer- 
entially, a massive companion star (with M > 15 Mq), given 
the dense environment in which it is expected to form and 
live, (iii) experiences a phase of wind fed or RLOF accretion 
that terminates when the donor star ends its life, and that 
(iv) the donor star is sufficiently massive to explode leaving 
behind a radio pulsar, as a remnant. 

Hypothesis (iii) implies that an X-ray active phase pro- 
ceeds the formation of the IMBH-PSR binary. Thus, IMBH- 
PSR progenitors can be identified as X-ray sources. The X- 
ray active phase is found to be characterized by luminosi- 
ties that spread over a broad range, depending on the mass 
of the companion star and the initial orbital separations. 
Two possibilities may occur: (a) If the star on the ZAMS 



is underfilling its Roche lobe, WFA luminosities in excess of 
iuLX ~ 3 x 10 39 erg s _1 are found only in presence of donors 
as massive as 30 Mq. With lighter stars, WFA luminosities 
fall typically in the range between 10 37 — 10 39 ergs -1 , with 
the highest values coming from transfer rates magnified by 
the focussing action of gravitational field of the heavy IMBH 
on wind particles. At separations upto several AU, typical 
WFA luminosities fall short below 10 38 erg s _1 so that, based 
on luminosity arguments, it would be impossible to distin- 
guish an accreting IMBH from a stellar-mass black hole, 
under WFA. WFA can last as long as the entire time spent 
by the star on the main sequence, tms ~ 1 — 10 Myrs. Ex- 
pansion of the massive star during nuclear evolution leads 
inevitably to a phase where the star fills its Roche lobe. This 
occurs when the binary separation a £ 10 AU at the time the 
star is transiting the Hertzsprung gap and/or ascending the 
giant branch. At this moment, one expects a dramatic in- 
crease of the mass transfer rate (see the spikes in Fig. 3 and 
4, and also Portegies Zwart et al. 2004b). This is a short 
lived phase of a few 10 4 — 10 5 yrs that accounts only for 
~ 5% of tms, in which the binary can emit luminosities in 
excess of 10 40 ergs -1 , characteristic of those bright ULXs in 
which we can be confident to find an IMBH. Thus, the life- 
time of a ULX, progenitor of an IMBH-PSR system, txjlx 
can be as short as a few 10 4 yrs. 

(b) If, instead, the massive star overfills its Roche lobe near 
the zero age, contact is maintained stably in the binary, and 
mean transfer rates as large as 10 -7 — 10 -5 Mq yr -1 are 
found. These rates are large enough to guarantee the stabil- 
ity of the disc against thermal-ionization perturbations. For 
donors more massive than ~ 15 M© the mass transfer rate 
is higher than the Dubus' critical value so that we argue in 
favor of persistent emission in binaries with an IMBH and 
massive donors. The accretion luminosities inferred at an 
averaged efficiency of 0.1, all exceed 10 39 ergs -1 , but only 
stars more massive than 20 Mq can guarantee luminosities 
above lO^ergs" 1 during the main sequence. RLOF with 
massive donors (> 20 M©) at the ULX level lasts for a time 
tulx of few million years, while lighter donors are found to 
brighten at the level of a few 10 39 ergs -1 for a time of > 10 
Myrs. 

Concerning the ultimate fate of the donor star we have 
seen that RLOF on the ZAMS causes significant mass losses, 
so, one can estimate that only stars heavier than 15 Mq and 
lighter than 30 Mq may end their lives as neutron stars and 
thus turn on as radio pulsars. Lighter donors end proba- 
bly their life as white dwarfs. These can still produce ULX 
activity under transient conditions. 

With hypothesis (iv), the formation of an IMBH-PSR 
system requires that the supernova explosion of the com- 
panion to the IMBH does not end with the disruption of 
the binary. Symmetric mass loss would not unbind the sys- 
tem given the large inertia of the IMBH, but the occurrence 
of an asymmetric natal kick can place the neutron star in 
an unbound orbit. Considering the pre-explosion orbital pa- 
rameters from our binary evolution models, we find that the 
probability of survival is very high, around 99%. The IMBH- 
PSR binary that survives disruption would have preferential 
separations between 1 and 10 AU, and a mild eccentricity. 
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5.2 Perspectives of detection of the radiopulsar 

In our study we have shown that binaries hosting an IMBH 
and a donor star, massive enough to leave a pulsar as a relic, 
likely experience an X-ray phase characterized by accretion 
luminosities ranging from well below the Eddington lumi- 
nosity of a one solar mass star, up to values characterizing 
the bright tail of the ULX window. The life time tulx over 
which a ULX phase is observed falls in the range between 
0.05 to 10 Myrs whereas rx ^ 10 Myrs is the characteris- 
tic lifetime of an accreting IMBH as "normal" X-ray source. 
While the formation of a binary IMBH-PSR system does not 
necessarely require a transit through a ULX phase, a ULX 
phase can inversly highlight a system that will evolve into 
an IMBH-PSR binary. From the knowledge of the distribu- 
tion function P(a) of the initial orbital separations we can 
calculate the ratio 5ft > 1 between the number of binaries 
hosting an IMBH-massive star system ending as a IMBH- 
PSR, to the number of binaries ATulx which experience a 
bright ULX phase before becoming an IMBH-PSR. 5ft is: 



5ft = 



/° m " P(a)da 
/" ULX P(a)da 



(15) 



where r t is the tidal distruption radius of the IMBH, a max is 
the orbital separation at which an X-ray phase (of arbitrary 
low intensity) occurs ending with the formation of a pulsar 
bound to the IMBH, and oulx is the initial orbital separa- 
tion at which a ULX phase sets in during binary evolution. If 
we suppose P(a) is uniform, and take a max ~ 200 AU as the 
limiting distance above which a natal kick would unbind the 
system (for a chance probability of 50%), and auLX — 1AU 
(the distance for a RLOF phase during the main sequence) , 
then we obtain an approximate estimate of 5ft ~ 200. This 
shows that the IMBH-PSR binaries whose progenitor ex- 
perienced a ULX phase may well represent only the tip of 
the iceberg of the whole population of IMBH-PSR binaries. 
However, we note that the estimate of 5ft is strongly affected 
by the uncertainity of the maximum orbital separation al- 
lowed by an exchange interaction (which corresponds to the 
cutoff of the distribution P(a)), and also by the real form 
of the distribution P(a). Detailed stellar dynamical simula- 
tions are required to improve the knowledge of the crucial 
factor 5ft. 

Provided 5ft is known, we can estimate the num- 
ber of IMBH-PSR binaries, from the simple formula 
(tpsr/tx)5R/A'ulx where tpsr ~ 10 Myr is the lifetime of 
a radiopulsar, A^ulx the observed number of bright ULX 
sources and / < 1 is a fudge factor to compute A^lx = 
/A^ulx from the known number of ULXs. In order to cal- 
culate the number of potentially observable radiopulsars 
ATimbh-psr orbiting an IMBH, one must also account for 
the strong anisotropy in pulsar emission. Introducing the 
pulsar beaming factor b ~ 10 (Lorimer 2001), and scaling 
the others quantities to the reference values, it turns out 



A/imbh-psr —20/ 



Tpsr 
10 Myrs 



/ 10 Myrs \ 5R 10 

V 7X / 200T 



Nil 



(16) 



In section §4 we have shown that Doppler distortion of the 
radio signal should not have hampered the detection of a 
putative IMBH-PSR system by the Parkes Multibeam sur- 
vey of the Galactic Disc (Manchester et al. 2001). However, 
other selection effects influence the detectability of a given 



radio pulsar during a survey, the most relevant being the 
pulsar intrinsic luminosity and the survey sensitivity limit 
for a source at the given position in the Galaxy, in turn 
depending on the pulsar spin period, duty cycle, dispersion 
measure, and amount of scattering smearing. A population 
synthesis analysis combined with a simulation, for the given 
survey, is necessary (e.g. Lorimer et al. 1993) for assessing 
the ratio x ^ 1 between the number of expected detections 
and the number of potentially observable radio pulsars (i.e. 
those having the radio beam sweeping the line of sight to 
the Earth). Putative pulsar companions to IMBH and typ- 
ical isolated pulsars are populations of non-recycled pulsars 
born and still residing in the Galactic disc. Assuming, in 
first approximation, that the two populations both share 
the galactic distribution and the same intrinsic properties 
(such as the luminosity function), we can apply to the pul- 
sars orbiting an IMBH the "average" value of \ properly 
calculated for the Parkes Multibeam Survey (Manchester et 
al. 2001). In particular, \ is in the range 0.050 — 0.075 for 
pulsars with luminosity greater than 1 mjy kpc 2 (Vranese- 
vic et al. 2004), and is ~ 3 times smaller, when including 
fainter sources having minimum luminosity ~ 0.2 mjy kpc 2 . 

Correcting for this factor, and adopting the reference 
values of equation (16), we expect that the Parkes Multi- 
beam survey should have detected 



A^obs = XNlMBH-PSR — /A?ULX 



(17) 



pulsars orbiting an IMBH. This estimate implies that the 
detection of a radio pulsar orbiting an IMBH in the Galaxy 
is subjected to the same small number statistics of the ULX 
population (whose detection in the X-ray band is free from 
biases), even for the most favorable case / = 1. In particular 
it is compatible with the current lack of observed pulsars of 
this type. More sensitive surveys of the Galactic disc (such 
as the ongoing P- Alpha survey at Arecibo) or targeted sur- 
veys (searching deeply for pulsars in the directions of young 
star clusters) may results in a improved value of x, enhanc- 
ing the possibility of unveiling a pulsar orbiting an IMBH 
within few years. However another possibility is that the 
number of ULXs hosting an IMBH is not so high as sup- 
posed, since as previously discussed, these sources may be 
explained by simple stellar mass black holes, leaving only 
a few number of ULXs with a true IMBH. In this case we 
must rescale the value of the equation (16) to the unknown 



number N rj 



< Nu 



of ULXs with a real IMBH, dimin- 



ishing our chance of a successfull detection. Next generation 
radiotelescopes will be able to go well beyond the Magel- 
lanic Clouds and eventually reach all the galaxies in the 
Local Group (e.g. M31 with LOFAR, Rottgering 2003) or 
nearby galaxies (e.g. M82 with SKA, Cordes et al. 2004). 
The much larger volume of space explored will enclose some 
known bright ULX 1 . If these future instruments will be able 
to sample a sizeable fraction of the total pulsar population 
in nearby external galaxies, they will greatly increase the 
chance of discovering objects belonging to the very intrigu- 
ing class of binaries presented in this paper. 



1 The nearest known bright ULXs are at a distance of ~ 2— 3Mpc, 
whereas the ULX reported in the spiral galaxy M33 seems to be 
an ordinary X-ray binary in a high state rather than a true ULX 
(Foschini et al. 2004). 
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